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Entanglement dynamics of two independent qubits 
in environments with and without memory 



o 
o 

> 

o 

G\ 
(N 

Oh' 



> 



o 



X 



B. Bellomo, R. Lo Franco, and G. Compagno 
Dipartimento di Scienze Fisiche ed Astronomiche, 
Umversita di Palermo, via Archirafi 36, 90123 Palermo, Italy 
(Dated: February 2, 2008) 

A procedure to obtain the dynamics of independent qudits (d-level systems) each interacting 
with its own reservoir, for any arbitrary initial state, is presented. This is then apphed to study 
the dynamics of the entanglement of two qubits, initially in an extended Werner-like mixed state 
with each of them in a zero temperature non-Markovian environment. The dependence of the 
entanglement dynamics on the purity and degree of entanglement of the initial states and on the 
amount of non-Markovianity is also given. This extends the previous work about non-Markovian 
effects on the two-qubit entanglement dynamics for initial Bell-like states [B. Bellomo et al., Phys. 
Rev. Lett. 99, 160502 (2007)]. The effect of temperature on the two-qubit entanglement dynamics 
in a Markovian environment is finally obtained. 
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I. INTRODUCTION 

One of the most striking characteristics of quantum 
mechanics is represented by the entanglement. It plays 
a role in fundamental aspects, like non- locality [ij, 
and applicative ones, like quantum information process- 
ing 3 . It is therefore of interest the experimentalgener- 
ation and manipulation of entangled systems [1, 0| 
and the theoretical study of the entanglement evolution. 
It is an important subject the analysis of entanglement 
decay and its relation with decoherence induced by the 
unavoidable interaction of the systems of interest with 
the environment. 

The study of the entanglement evolution of two qubits, 
initially in a Bell state and exposed to local zero temper- 
ature and Markovian noisy environments, has led to the 
surprising result that it may markedly differ from the 
single qubit coherence evolution. In particular, in front 
of an exponential decrease of the single qubit coherence, 
the entanglement may instead completely disappear at a 
finite time [1, [1] , a phenomenon known as " entanglement 
sudden death" (ESD) which has been proven to occur in 
a quantum optics experiment [To| . ESD appears to put a 
limitation to the time when entanglement could usefully 
exploited. Its intrinsic and practical interest has led to 
several investigations in Markovian environments at ei- 
ther zero or finite temperature and for initially pure and 
mixed states [HI [H [ll H [H [H, [13, [3 . The disen- 
tanglement phenomenon being linked to the creation of 
body-environment correlations makes interesting to ana- 
lyze its behavior in a non-Markovian system where mem- 
ory effects are relevant 

Recently, a procedure has been developed that allows 
to obtain the complete dynamics of the two-qubit system 
by knowing the dynamics of each qubit interacting with 
its environment, which is applicable to any initial condi- 
tion. This has permitted to show, for two noninteracting 
qubits initially described by pure Bell-like states and in 
non-Markovian environment at zero temperature, that 



entanglement may reappear after some period of time of 
complete disappearance [2^ . This appears relevant from 
a fundamental point of view and because also may permit 
to extend the time of usefulness of entanglement. 



The aim of this paper is to generalize the results ob- 
tained in Ref. ^] first giving a formal extension of the 
procedure described above from two qubits to N indepen- 
dent qudits (d- level systems), each interacting with its 
own reservoir. Then, we shall apply it to determine the 
entanglement dynamics of two identical qubits, initially 
described by a generalization of mixed Werner- like states, 
in non-Markovian environment at zero temperature, ex- 
tending to a wider class of initial conditions the previous 
study. This will permit to find how entanglement dynam- 
ics and its revivals are related to physical parameters like 
the purity and the degree of entanglement of the initial 
states or to the amount of non-Markovianity of the sys- 
tem. Then, we shall exploit our procedure to extend in a 
simple way the previous analysis of two-qubit entangle- 
ment in a Markovian environment at finite temperature 
to the initial extended Werner-like states. This will per- 
mit to link the occurrence time of ESD to the purity or 
the degree of entanglement of the initial states. 

The paper is organized as follows. In Section |TT] we 
propose a general procedure to solve the dynamics of 
N independent qudits. In Section |TTT] we give a stan- 
dard Hamiltonian model describing the dynamics of sin- 
gle qubit-reservoir interaction, and we also determine the 
two-qubit density matrix evolution. In Section ITVl we in- 
troduce our initial states and recall the definition of con- 
currence to quantify the two-qubit entanglement. In Sec- 
tion |V] we study the two-qubit entanglement dynamics in 
non-Markovian environments at zero temperature while 
in Scction lVTl we analyze the effects of the temperature on 
the entanglement dynamics in Markovian environments. 
In Section [VIII we summarize our conclusions. 
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II. PROCEDURE 

We consider a system formed by N non-interacting 
parts S = 1,2, . . . , N, each part consisting of a d level 
system (qudit) S = 1,2, ... ,N, locally interacting re- 
spectively with a reservoir Rs — Ri, R2, . . . , Rn- Each 
qudit S and the corresponding reservoir Rg are initially 
considered independent. The evolution of the reduced 
density matrix of the single qudit S is given by 

p'{t)=Trns {uHt)pHO)^p''HO)U'Ht)}: (1) 

where the trace is over the reservoir Rs degrees of free- 
dom and U^{t) is the time evolution operator for the 
part S. By using the spectral decomposition of the den- 
sity matrix p^^ (0) of the reservoir 

p«-(0) = ^A.,|^.,)(¥..,|, (2) 

as 

the reduced density matrix p^{t) can be expressed in the 
so called operator-sum representation as [3| 

pHt) = E w!sPsit)p'i(^)wl!,At), (3) 

where the operators W^_^ (t) are given by 

W^sPs = \/^(VPs\uHt)\Vc.s)- (4) 

The assumption of independent parts implies that the 
time evolution operator U (t) of the complete system 1 -|- 
2 + ... + N factorizes as U{t) = U^{t) (g) C/^ (i) (g) . . . (g) 
U^{t). It follows that the operator-sum representation 
of the reduced density matrix for the A'^-qudit system 
{1 + 2 + ... + N) reads like 

^p^im^f3AtwiX(t) ■ ■ ■ w^^jt). (5) 

Given the basis {\is),is = ^S^'^s, ■ ■ ■ ids} for the qudit 
S, inserting the unity operators Is — J2 i^sl between 
operators and density matrix in Eq. ([3]), it follows that 
the dynamics of each qudit has the form 

PMiW = E<5w/^Mi(o)' 

ph'S') = E4:twpLi(o) 
12V2 

pL'Jt) - E ^'t^pti^'s^)- (6) 

Adopting the same procedure for the A^-qudit reduced 
density matrix p{t) given in Eq. the dynamics of the 



iV-qudit system is correspondingly expressed by 

h^l ^2^2 ^Af^jv 

where is = Is, 2s,..., ds with S ^ 1,...,N. Eqs. ® 
and ([7|) clearly show that the dynamics of A^-qudit den- 
sity matrix elements follows by knowing the ones of each 
single qudit, and that it is applicable to any initial con- 
dition of the total system. 

We wish to point out that the procedure developed 
above, which permits to obtain the dynamics of N in- 
dependent qudits each locally interacting with its own 
reservoir, can be also applied to the case where the reser- 
voir is unique but the qudits are placed at distances larger 
than the spatial correlation length of the reservoir. An- 
other point we wish to stress is the practical aspect of 
our procedure. In fact, it allows to obtain the dynamics 
of N qudits, provided that the dynamics of each qubit 
is known, by purely algebraic way and independently 
from the initial conditions. This considerably simplify 
the work typically required by writing and solving differ- 
ential equations for the density matrix elements of the 
total system [2^ or by using other approaches whose 
difficulty increases with the number of parts, levels or 
environmental noise complexity of the total system [l^ . 



III. MODEL 

We now shall apply the above procedure to the case of 
a system composed by two parts, each one consisting of a 
two-level system (qubit) interacting with a reservoir. We 
shall take the single part "qubit -f reservoir" described by 
the Hamiltonian 

H = hwoa+a- + E ^fi^kblbk + {gkcr+bk + gla^V^ , 

k 

(8) 

where loq is the transition frequency and a± are the sys- 
tem raising and lowering operators of the qubit, while 
bj, , bk are the creation and annihilation operators and gk 
the coupling constant of the mode k with frequency ujk- 
The Hamiltonian of Eq. ([8]) is used to describe various 
systems such as a qubit formed by an exciton in a po- 
tential well environment. However to fix our ideas we 
shall take it to represent a qubit formed by the excited 
and ground electronic state of a two-level atom interact- 
ing with the reservoir formed by the quantized modes of 
a high-Q cavity. The dynamics of an atom interacting 
with a quantum electromagnetic field under the dipole, 
rotating-wave, and two-level approximation is in fact de- 
scribed by the Hamiltonian of Eq. ([H]). 

Systems described by this Hamiltonian have been stud- 
ied under different conditions, as weak coupling and finite 
temperature in both Markovian [2^ and non-Markovian 
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[27| approximation, and strong coupling at zero temper- 
ature [H, [2^. In all of these cases, as we shall report in 
Sections |V] and IVH the single-qubit density matrix evo- 
lution can be cast in the following form: 



pfiit) pUt) 

Pilit) Ploit) 



(9) 



where 



ploit) = (l-uf)pfi(0)-f (l-i;f)p^o(0), 

pfoW = pEiit) = ^fpfo(O), (10) 

uf^vf^zf being functions of the time t. Now we are 
ready to use, following the procedure described before, 
the evolution of the reduced density matrix elements for 
the single qubit to construct the reduced density matrix p 
for the two-qubit system. The two qubits, labeled by 5* = 
A, B, may be in general in different environments so that 
their evolution is characterized by different values of the 
functions wf,wf,zf. In the standard basis B — {|1) = 
111), |2) = |10), |3) = |01), |4) = |00)}, using Eqs. ©, © 
and (jlOp. we obtain for the diagonal elements 

pil(t) = ufufpil{Q)+utvfp22{Q) 

-f v^ufp^m + vtvfp^M, 

= Wf(l-uf)pil(0)-|-uf(l-«f)p22(0) 
+ vtil - uf )p33(0) + vt{l - vf)p44{0), 

P33{t) - {l-uf)ufpn{0) + {l-ut)vfp22{0) 
+ (1 - vf)ufp,^{0) + (1 - vf)vfp^,{0), 
PM{t) = {l-ut){l-uf)pn{0) + (l-ut) 

X (1 - Vf)p22{0) + (1 - vf){l «f )P33(0) 



+ {l-vf){l-vf)p,,{0), 
and for the non-diagonal elements 



(11) 



Pl2{t) 

Pisit) 

Pli{t) 
P2i{t) 
Pii{t) 



utzfpi2{0)+vfzfp3m 
zfufpi3{0) + zfvfp2m 

zfzfp.^iO), p23{t) 



P23(0), 



zfil ~ )pi3(0) zt{\ - vf)p2M. 

(1 - «f )zf Pi2(0) (1 - ^;^^)zf P34(0),(12) 
with Pij{t) ~ Pji{t), p{t) being a hermitian matrix. We 



point out that Eqs. (fTT|) and (fT^ . being valid for any 
form of the single-qubit evolution functions uf , zf , 
are more general of the corresponding ones given in 
Ref. [13] , where the reduced density matrix elements have 
been given for the particular case of non-Markovian en- 
vironment at zero temperature, and they also allow to 
study the two-qubit density matrix evolution for any ini- 
tial state. 



IV. INITIAL STATES AND CONCURRENCE 

The two-qubit entanglement dynamics has been pre- 
viously analyzed taking as initial states pure Bell-like 



states [23]. However, as said before, it appears of in- 
terest to obtain the entanglement dynamics of two inde- 
pendent qubits, each locally interacting with a reservoir, 
in the case of more general initial conditions and in par- 
ticular for Werner states II4I. The procedure we have 



developed is applicable to any two-qubit initial state and 
therefore also to Werner states. In particular, in this sec- 
tion we consider as two-qubit initial states an extension of 
the Werner states. These "extended" Werner-like (EWL) 
states are mixed states that may reduce to Werner mixed 
states [lO, HH, HI] or to Bell-like pure states. They shall 
allow us to put in evidence the effects of both the mixed- 
ness and the degree of entanglement of the initial states 
on the entanglement dynamics. 



A. Initial states 

It is known that Bell states and Werner mixed states 
are important in quantum information and computation 
processing [1, [s^. Werner- like mixed states are also 
shown to occur during an experimental generation of Bell 
states [13] • However, from Eqs. (fTTj) and it is evident 
that the two-qubit density matrix at time t is valid for 
arbitrary initial states and its structure crucially depends 
on their form. We use this property to study the two- 
qubit entanglement dynamics starting from EWL states 
defined as 



/}*(0) = r]<i>)(cf>] 



1 - r 



/5*(0) = r]*)(vI/] + 1-^/4. 



(13) 



In Eq. (I13p . r indicates the purity of the initial states, li 
is the 4x4 identity matrix and 

]$) =a]01)+/3]10), ]vl') = a]00) + /3]ll), (14) 

are the Bell-like states with a real, fi = ]/3]e"^ and 
Q;2-h]/3]2 = 1. The states defined by Ea, (HSl) reduce to the 
well-known Werner-like states [3ll . l32| when their pure 
part becomes a Bell state, that is for a = ±/3 = 1/^/2. 
The degree of entanglement of the Werner-like states can- 
not be increased by any unitary transformation, so that 
they correspond to a maximum degree of entanglement 
for a given linear entropy. For r — Q the Werner-like 
states become totally mixed states, while for r = 1 they 
are the well-known Bell states. On the other hand, fixing 
r — \, the EWL states of Eq. p3|) reduce respectively to 
the Bell-like pure states \^),\^) of Eq. dH]). The EWL 
states /5*(0),/5*(0) are therefore more general than the 
Werner-like states since they contain the possibility that 
their pure part is not maximally entangled. They, as 
said before, thus permit to study the effect of both the 
initial states mixedness and the degree of entanglement 
of their pure part on the entanglement dynamics. The 
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initial density matrix elements are for the state 
l-r * 1 



pfi(O) 



Pf3(0) 
Pf4(0) 



4 

1 2 

+ a r. 



Pl2(0) 



4 

0, Pf3(0) 
and for the state 

l-r _ 



4 

l-r 



pKo) = 

(0) = a/3r, 



1 - r 



+ 
1 - 



(15) 



^4*4(0) = 

Pli0) = 



The structure of the density matrix elements of Eqs. p5)) 
and show that the EWLSs belong to the class of 
the "X" states, that is to those states having non-zero 
elements only along the main diagonal and anti-diagonal. 
The general structure of an "X" density matrix is thus 



/5(0) 



/ Pii(O) pi4(0) 

P22(0) /923(0) 

P23(0)* P33(0) 

Vpi4(0)* P44(0) 



(17) 



Bell states and Werner states belong to this class of states 
and we also point out that a "X" structure density matrix 
may also arise in a wide variety of physical situations ^34 . 
IHjIUjIs^I- a remarkable aspect of the "X" states is that, 
under the single qubit evolutions of Eq. (jlOp determined 
by the Hamiltonian of Eq. ([8]) , the initial "X" structure 
is maintained during the evolution, as easily seen from 
Eqs. pT|) and Moreover, it is know that, for such a 
"X" density matrix, the expression of the concurrence in 
terms of the density matrix elements takes a particular 
simple form [T^ . 



B. Concurrence 

In order to describe the entanglement dynamics of the 
bipartite system, we use Wootters concurrence |38[ . This 
is obtained from the density matrix p of the two-qubit 
system ^ -I- _B as 



{0, VAi - VA2 - VA3 - VA4}, (18) 



where the quantities are the eigenvalues of the ma- 
trix C = piPy ® cry)p*{ay (g) ay), arranged in decreasing 
order. Here cr^ are the well-known Pauli matrices and 
p* denotes the complex conjugation of p in the standard 
basis B. Cp{t) varies from for a disentangled state to 
1 for a maximally entangled state. In our case, Cp(t) as 
a function of the two-qubit density matrix elements of 
Eqs. (fTTj) and ^2]) depends on the initial condition and 
on the dynamics as expressed by the functions wf , , . 

The concurrence at the time t for an initial "X" state 
defined in Eq. ([T7|) can be easily computed by exploiting 



the fact that, under our dynamical conditions, the "X" 
structure is maintained during the two-qubit evolution. 
It results to be given by [l3| 



2ma.Ji{0,Ki{t),K2{t)}, 



where 



Ki{t) = |p23(i)| - v/piiWp44(i) 

/^2(<) = |Pl4(i)| - V/P22(i)p33(i) 



(19) 



(20) 



The EWL states of Eq. (|13p have the same initial value 
of the concurrences, given by 

(16) C*(0) = C*(0) = 2max{0,(a|/3Kl/4)r-l/4}, (21) 



from which one finds that there is initially entanglement 
when r > r* — (l-|-4a|/3|)~^, thus entanglement exists for 
these states only if r > 1/3. In particular, the Bell-like 
pure states |$), |*), obtained by p*(0),p*(0) for r = 1, 
have the same degree of entanglement (concurrence) reg- 
ulated by a according to C$(0) = C*(0) = 2a\J\ — o? . 
Using the expression of the concurrence for "X" states 
given in Eqs. (fTOjl and PU)) . it is readily seen that C^{t) 
for the EWL states of Eq. (fT3|) is respectively 



Cf{t) = 2max{0, Ki{t)}, (t) = 2max{0, K2{t)}. 

(22) 

We shall use these equations in order to study the two- 
qubit entanglement dynamics. In particular, from now 
on we shall limit our analysis to the case of two identi- 
cal qubits locally interacting with identical environments. 
This implies that the functions uf , zf {S = A,B), in- 



troduced in Eq. (|T0)) . are the same for both qubits and 
we shall then indicate them simply by ut,vt, zt- 



V. NON-MARKOVIAN ENTANGLEMENT 
DYNAMICS 

In this section, we shall analyze the two-qubit entan- 
glement dynamics for our EWLSs of Eq. in the non- 
Markovian case at zero temperature. This constitutes an 
extension of our previous study of the non-Markovian ef- 
fects on the two-qubit entanglement dynamics for the ini- 
tial Bell-like pure states of Eq. [24] . This extension 
shall permit to point out how the two-qubit entangle- 
ment dynamics is affected by the presence of mixedness 
in the initial states regulated by the purity parameter r 
and the degree of entanglement represented by a. On 
the basis of our procedure described in Sec. |TT1 we shall 
obtain the reduced two-qubit dynamics by the knowledge 
of the reduced single qubit dynamics. 

Single qubit dynamics induced by the Hamiltonian of 
Eq. ([5]) in the non-Markovian limit has been treated in 
several works [13, [H,!!^. In the following, we exploit the 
results there obtained about this dynamics under weak 
coupling at finite (low) temperature and strong coupling 
at zero temperature. 
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A. Weak coupling (low T) 

The exact single qubit dynamics induced by Hamil- 
tonian of Eq. ([5]), with the field initially at finite tem- 
perature, has been recently studied [271], and a formal 
but non transparent, because of its complexity, solution 
for the density matrix elements evolution has been ob- 
tained. However, at low temperature (g-'^'^o/fceT ^ j^-j 
and weak coupling (F/wq ^ 1, T being the zero tem- 
perature spontaneous emission rate), a simple explicit 
solution has been found, that can be put in the form of 
Eq. (Uni) with 



Ut = 1- 



Vt 



(l-e-i*)(l-e" 



(1 



■^0 



1 - e 
1 — e ''B^ 



^^0 _p+ 
fc „ T 



Zt = 



1 



-rt/2-iuJot 



rt 



(23) 



This solution presents some differences with respect to 
the corresponding low temperature Markovian dynam- 
ics, for example the upper level relaxation time results 
longer than the Markovian one [l^]. However, it is pos- 
sible to show that the evolution of the concurrence for 
the two initial states /5*(0), jO*(0) of Eq. (fT5|) is not qual- 
itatively much different from the Markovian one that we 
shall treat in Sec lVIl In the low temperature regime, a 
strong qubit-reservoir coupling condition seems to be re- 
quired in order to have evident non-Markovian effects in 
the entanglement dynamics of a two-qubit system. To 
this aim, in the following we shall then use the results 
for the non-Markovian strong coupling regime obtained 
at zero temperature j29j . 



B. Strong coupling (T = 0) 

Following our previous work [23 | , here we use the exact 
solutions for the single qubit dynamics at zero tempera- 
ture [2^ , [2^ with the aim to extend the analysis of two- 
qubit entanglement dynamics to the more general initial 
conditions assigned by the EWLSs of Eq. p3)) . 

Because the Hamiltonian of Eq. ([5]) represents a model 
for the damping of an atom in a cavity, in the case of a 
single excitation in the atom-cavity system, the effective 
spectral density J{uj) is taken as the spectral distribution 
of an electromagnetic field inside an imperfect cavity sup- 
porting the mode luq, resulting from the combination of 
the reservoir spectrum and the system-reservoir coupling, 
with r related to the microscopic system-reservoir cou- 
pling constant. This spectral density has the form [2^ 



Jico) 



1 



(24) 



where A, defining the spectral width of the coupling, is 
connected to the reservoir correlation time tb by the rela- 
tion tb ~ A^^. In fact, it can be shown that the reservoir 
correlation function, corresponding to the spectral den- 
sity of Eq. (|24p . has an exponential form with A as decay 
rate. On the other hand the parameter F can be shown 
to be related to the decay of the excited state of the atom 
in the Markovian limit of flat spectrum. The relaxation 
time scale tr over which the state of the system changes 
is then related to F by th « F~^. The exact solution 
for the single-qubit density matrix at zero temperature 
with the spectral density given by Eq. has the form 
of Eq. (Uni) with [H,!!^ 



Ut 


= e-^' 


Vt 


= 0, 


Zt 





dt 

cosl- 



A . fdt 
- sm — 
d \ 2 



(25) 



where d — y/2TX — A^. A weak and a strong coupling 
regime can be distinguished for the single qubit dynam- 
ics. The weak regime corresponds to the case A/F > 2, 
that is Tn > 2tb- In this regime the relaxation time is 
greater than the reservoir correlation time and the dy- 
namics is essentially Markovian with an exponential de- 
cay controlled by F. In the strong coupling regime, that 
is for A/F < 2, or t^j < 2tb, the reservoir correlation 
time is greater or of the same order of the relaxation 
time and non-Markovian effects become relevant. There- 
fore, we shall study the two-qubit entanglement dynam- 
ics in this regime. Substituting the expressions for ut, 
Vt and Zt given by Eq. ([25]) in Eqs. ([TT]) and ([12]), the 
two-qubit density matrix at time t for both the initial 
states of Eq. is obtained. Using for the concurrence 
Eq. (j22p . we then determine its dependence on time i, 
purity parameter r, probability amplitude a and a di- 
mensionless parameter A/F that indicates the coupling 
regime. The explicit expressions of concurrence in terms 
of these parameters are, although simple, rather long and 
not particularly enlightening. Therefore, we shall plot 
the concurrence as a function of both Tt and one of the 
dimensionless parameters (r, a^, A/F), after fixing the re- 
maining ones. 

The case for r — 1, that is when the initial EWL states 
of Eq. (Uni) reduce to the pure Bell- like states of Eq. pi)) , 
has been previously treated [2j]. The concurrence time 
evolution in terms of Ft and has been analyzed, for 
the fixed value A/F = 0.1. This choice of A/F corre- 
sponds to an experimentally feasible value in the cavity 
QED context. In fact, cavity QED experimental configu- 
rations have been recently realized using Rydberg atoms 
with lifetimes Tat ~ 30ms, inside Fabry-Perot cavities 
with quality factors Q « 4.2 x 10^" giving a cavity life- 
time Tcav ~ 130ms 40], these values corresponding to 
A/F « 0.1. In particular, it has been shown that C$(t), 
related to the initial Bell-like state /5*(0) = 1<I')(<I'1, pe- 
riodically vanishes according to the zeros of the function 
Ut of Eq. (j25p with a damping of its revival amplitude 



FIG. 1: (Non-Markovian case; r = 1, = 1/3, T = 0, A/F = 
0.01) Comparison between concurrences Ci^{t) (dashed line) and 
C\ii {t) (solid line) as a function of the dimensionless quantity Tt. 



FIG. 2: (Non-Markovian case; T = 0, a'^ = |/3p = 1/2, A/F = 0.1) 
Concurrence Cp(t) (left figure) and (t) (right figure) for Werner- 
like states as a function of the dimensionless quantities Tt and r. 



and ESD does not occur. The concurrence C>i,(i), re- 
lated to the initial state /5*(0) = has been shown 
to have a similar behavior for values of > 1/2, while 
for < 1/2 two ranges of the parameter may be dis- 
tinguished, one where ESD occurs and another where 
revival of entanglement appears after periods of times 
(dark periods) when entanglement has completely dis- 
appeared. As an example of this behavior, we plot in 
Fig.[T]the concurrences C$(t) and C>i,(t) for — 1/3 and 
X/T — 0.01 which corresponds to a strong non-Markovian 
coupling. This choice of X/T is taken to make the revival 
phenomenon evident. Fig. [T] shows that the remarkable 
phenomenon of entanglement revivals occur, for the state 
/5*(0) = after finite dark periods. Moreover, we 
wish to point out that, for initial Bell states, the concur- 
rences at the time t take the simple form 



C^{t) = ut, C^it)^ui, (r = 1, a ==1/31 = 1/ V2 



(26) 



with the zeros of the function ut of Eq. ((25l) given by 



tn = —[mr — arctan((i/A)]. 
a 



(27) 



The distance At between two consecutive zeros of Ut and 
then of the concurrences of Eq. is constant, depend- 
ing only on the ratio A/F as 



At 



27r/F 



v/A/r(2 - X/T) 



(28) 



Because our procedure can be easily applied to any 
initial state, we extend the analysis of the two-qubit en- 
tanglement dynamics for non-Markovian environments 
from pure Bell-like states to the interesting case of initial 
mixed states 12] of the EWL form given by Eq. (fT5)) . 
Therefore, the new aspects we wish to examine here 
correspond to the case when the purity r of the initial 
EWL states of Eq. (|13p is less than one. In particular, 
the concurrences C* (t) and C* (t) are plotted in Fig. [2] 
in terms of the dimensionless quantities Tt and r, for 
X/T = 0.1 and a = l/V^. This choice of a corresponds 
to take Werner-like states as initial states, as described 
in Sec. IIV Al Different choices of a values do not give 




FIG. 3: (Non-Markovian case; T = 0, r = 1, = 1/3) Concur- 
rence Cp{t) (left figure) and (t) (right figure) as a function of 
the dimensionless quantities Ft and A/F. 



entanglement dynamics qualitatively different from the 
case treated here. From Fig.[2]we see that ESD occurs for 
both the initial Werner- like states /5*(0) and /5*(0) and 
sufficiently small values of r. This is at variance with the 
behavior of initial Bell-like states, where ESD occurred 
only for the state /5*(0) — |\E')(^'|. For intermediate val- 
ues of purity r, entanglement revivals after finite dark pe- 
riods occur for both initial Werner-like states. This again 
differs from previous results [24] where revivals and dark 
periods of entanglement occurred only for the Bell-like 
state ]\1/). In the region of r large (r ^ 1) the time be- 
havior of the concurrences approaches the one of initial 
Bell states. The reason of the decrease of the concur- 
rences value at i = with the decrease of the purity r is 
linked to the corresponding increase of the initial states 
mixedness. On the other hand, the period of time of en- 
tanglement revivals increases with the increase of r and 
thus with the increase of the initial states purity. 

Another aspect of interest is how the entanglement dy- 
namics is influenced by the values of A/F, which regu- 
lates the degree of non-Markovian characteristics of the 
system. This is an aspect that has never been examined 
previously. To this purpose we plot, in Fig.[3l the concur- 
rences C* (t) and C* (t) as functions of the dimensionless 
quantities Tt and A/F, starting from pure Bell-like states, 
that is for r — 1 and taking a — l/VS- This choice of 
the parameters permits to evidence quite different time 
behaviors of the concurrences in terms of A/F for the two 
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initial Bell-like states p*(0) = |$)($| andp*(0) = |*)(*| 
of Eq. (fTil) . The reason of the choice a = 1/ \/3 is that 
when o? — 1/2, that is for initial maximally entangled 
pure states, the concurrences presents similar behaviors 
with the appearance of revivals of entanglement and dark 
periods but without the phenomenon of ESD. Instead, 
from Fig. [3] we see an evident difference in the entangle- 
ment dynamics of the two initial pure states. In partic- 
ular, starting from the pure state /5*(0) = |^')(^| with 
1$) = (|01) + \/2|10))/\/3, the left part of Fig. [3] shows 
neither ESD or dark periods of entanglement while the 
entity of the revival phenomenon increases by decreasing 
the coupling A/F, that is by enforcing the non-Markovian 
effects. On the other hand, starting from the pure state 
I*) = (|00) + %/2|ll))/V3, the right part of Fig. [3] shows 
ESD occurrence for high values of A/F, that is for small 
non-Markovian conditions, and revivals of entanglement 
after dark periods of time for sufhciently small values of 
A/F, where the non-Markovian effects are stronger. The 
same analysis for mixed initial states (r < 1) does not 
present a qualitatively different behavior of entanglement 
dynamics in terms of the amount of non-Markovianity 
regulated by the parameter A/F. 

These results summarize how the two-qubit entangle- 
ment dynamics in non-Markovian environment sensibly 
depends both on the initial conditions and the degree 
of non-Markovianity of the environments and show that 
different conditions may lead to very different time be- 
haviors of the entanglement. 



VI. MARKOVIAN ENTANGLEMENT 
DYNAMICS AT T 7^ 

We now shall analyze in detail the two-qubit entangle- 
ment dynamics in Markovian environments at finite tem- 
perature, starting from the initial EWL states defined in 
Eg. p^ and exploiting the well-known results obtained 
for the singlc-qubit reduced dynamics. In fact, on the 
basis of our procedure described in Sec. |TT1 we shall ob- 
tain the reduced two-qubit dynamics by the knowledge 
of the reduced single qubit dynamics. 

The single-qubit reduced dynamics associated to the 
Hamiltonian of Eq. ([5]) has been widely studied within 
Markov approximation for weak coupling. The single- 
qubit density matrix, p'^, satisfies the standard master 
equation in the Lindblad form [26| 



Ap^ F(n + 1) 



At 



(2(T_/5 (7+ — a+a^p — p a^aJ) 
^J(2a+p^a_ - a.a+p^ - p^a.a+), (29) 



where fi = l/(e — 1) is the thermal mean photon 
number of the resonant field mode at temperature T, F is 
the spontaneous emission rate and Fn the stimulated ab- 
sorption and emission rate. The solution of the previous 
master equation can be found by solving the system of 




FIG. 4: (Markovian case; = |/3p = 1/2, kBT/hwo = 0.1) 
Concurrence Cp(t) (left figure) and (t) (right figure) for Werner- 
fike states as a function of the dimensionless quantities Ft and the 
purity r. 



differential equations associated to each matrix element 
of p^ [2^ . The reduced density matrix elements can then 
be cast in the form of Eq. (|10[) , with the functions Ut , vt 
and zt given by 



Ut 



Vt 



Zt 



1 - e ''bt 
_e-rcoth(^)t-| 



-(rt/2) coth(. 



(30) 



Substitution of these expressions in Eqs. (fTTj) and (fT2|) 
allows to obtain the two-qubit density matrix at time t 
for both the initial states of Eq. (fTSll . The concurrence is 
then obtained by using Eqs. (PO)) and as a function 
of time t and the system parameters, that is purity r, 
probability amplitude a — -y/l — |/3p and temperature 
T. We recall that, at finite temperature, X states have 
been shown in general to alwa ys u ndergo sudden death 
of entanglement at finite time [25| . In order to describe 
how the ESD time does depend on the physical param- 
eters of the system. In the following we shall plot the 
concurrence as a function of Ft and one of the dimen- 
sionless parameters (r, a^, fc^T/fiiJo), for assigned values 
of the remaining ones. 



A. Low temperature 

The concurrences C*(t) and C*(t) of Eq. related 
respectively to the initial states /5*(0) and /5*(0), are 
plotted in Fig. [3] as functions of the dimensionless time 
Ft and purity r, for — \f3\'^ — 1/2 and for a given low 
temperature such that kBT/fkuo — 0.1. The choice of a 
corresponds, for the initial states, to the Werner-like ini- 
tial states, while the choice of low temperature rests on 
the fact that in this limit the two-level approximation for 
the atoms is more realistic. Fig. 2] shows that the ESD 
occurs for both states but at different times, the differ- 
ences depending on the values of purity r. In particular. 
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FIG. 5: (Markovian case; r = 1, kgT/hujo = 0.1) Concurrence 
Cp{t) (left figure) and {t) (right figure) for Bell-fike states as a 
function of the dimensionless quantities Tt and . 



FIG. 6: (Markovian case; r = 1, = 1 — |/3p = 1/3) Concurrence 
Cp(t) (left figure) and C^(t) (right figure) as a function of the 
dimensionless quantities Ft and kgT /Tloq. 



increasing the purity of the initial states, the ESD time 
is retarded for both states. Moreover, the ESD time for 
the initial state /5*(0) is longer than the one for /3*(0). 
This may be linked to the fact that the product of pop- 
ulations of the totally excited state (pii(t)) and ground 
state (p44(t)) of the two-qubit system can be shown to 
increase more slowly than the product of populations of 
the states where only one qubit is excited {p22{t), Pssit)), 
as seen from Eq. (1^ . 

For r = 1 the initial states reduce to Bell-like states, 
and we allow the degree of entanglement, represented by 
a^, to change again keeping ksT/hujo — 0.1. This differs 
from the T ^ cases studied previously, that were only 
restricted to initial Bell states [ll|. The concurrences 
Cf{t) and C*(t) are then given in Fig. O which shows 
that ESD occurs for both initial pure states of Eq. (fHl) 
independently on the initial degree of entanglement rep- 
resented by . This is at variance with the zero tem- 
perature behavior where ESD occurs only for the initial 
pure state p*(0) = |^')(*| and when a'^ <l/2^. The 
occurrence of ESD for both the initial pure states, when 
T 7^ and independently from the degree of entangle- 
ment, is due to the population with time of the qubit 
excited states. However, in the p*(0) case it is the popu- 
lation of only one of the qubit, represented by the terms 
P22{t) and P33{t) in Eq. (PO)) . that makes after a certain 
time K2{t) < in Eq. dH]). On the other hand, in the 
/5*(0) case it is the excitation of both qubits, represented 
by the term pii(t) in Eq. (I^Dl) . that makes after a given 
time Ki{t) < in Eq. The right part of Fig. [H 

shows for C* (t) an asymmetry in the times of ESD with 
respect to the maximally entangled case (a^ = 1/2), with 
shorter times for < 1/2. This is again attributable to 
an initially higher contribution of the upper states for 
a2 < 1/2 for the state p*(0) of Eq. (fn)). 

In order to show how temperature influences the evo- 
lution of concurrence, in Fig. [5] we plot the concurrences 
Cp (t) and C* (t) as functions of the dimensionless quan- 
tities Tt and kBT/hioo, fixing the purity (r = 1) and the 
initial degree of entanglement (a^ = 1/3). This choice of 
a allows to evidence a rather different behavior of con- 
currence for the two initial pure states /5*(0) = |^)($| 
and p*(0) — at small T. Moreover, for small 



values of kBT/fkuo, the ESD time is sensibly affected by 
slight temperature changes. Finally, for T = we also 
find that ESD occurs only for the pure entangled state 
p*(0) = |^')(*|, as previously obtained [ll|. 



B. High temperature 

The condition of low temperature (fc^T/ftwo <C 1) may 
result appropriate in some experimental contexts, as in 
the case of trapped atoms or ions [4l|,|42|. However, there 
may be experimental situations where the low tempera- 
ture condition cannot be applicable. This may be, for 
example, the case of two identical quantum dots each 
with energy difference between the first two levels given 
by E2 — El ^ ksT. Under the action of an external mag- 
netic field the ground state splits in two sublevels that 
may have an energy difference Hujq <C ksT. In this case 
one may apply, in the study of entanglement between 
the sublevels of two qubits, the two-level approximation 
also in what for these sublevels is the high temperature 
regime. 

In this high temperature limit (ksT /fno^ 3> 1), the 
concurrences of Eq. take a simple analytical form. 
Moreover, the concurrence C(t) takes the same form for 
both the two initial EWL states p*(0), p*(0) of Eq. (fT^ 
given by C{t) = max{0,^(i)} with 



(31) 



This expression is a generalization of the one obtained for 
pure maximally entangled states [l3| and takes into ac- 
count mixedness and degree of entanglement of the initial 
states. It leads to a dimensionless ESD time Fiosd whose 
expression in terms of temperature and initial degree of 
entanglement is 



Tt 



csd 




2a| 



(32) 



This expression shows that the dependence of t^sd on 
the initial purity and degree of entanglement is weak. 
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However, when (r - r*)/r* < 1 (r* = i+4a \ fj\ )' ^^^^ 
for small initial entanglement, one gets 

r?csd« 2^[(l + 4a|/3|)r-l] (33) 

and the ESD time has a linear dependence on the pu- 
rity r. When r substantially differs from r* , because of 
the weak dependence of tesd on r and a, we can limit to 
consider the case a = \l3\ — 1/V2 and r — 1 when the 
dimensionless ESD time of Eq. ([5^ assumes the partic- 
ularly simple form 

rtcsd« 0.882^ («' = |/3|' = l/2,r=l). (34) 

In typical experimental conditions, quantum dots are 
subjected to an external magnetic field _B ~ 1 - 10 T 
(43l.l43|. inducing a sublevel sphtting given by fkuo « fiB, 
where /x is the electron Bohr magneton, that corresponds 
to the dimensionless ESD time E^esd (3 x IQ-^ - 3)/T. 
For a temperature range T ^ (0.3 - 3)K we have 
tesd ^ 1/r. In this case, the limitation in the possible 
use of quantum gate operations with entangled quantum 
dots is determined by the ESD time instead that by the 
single quantum dot decoherence time. This lowers the 
number of operations M [4^ that the system can 
effectively perform. 

VII. CONCLUSIONS 

In this paper we have presented a procedure (Sec. |Tl| 
that allows to obtain the dynamics, represented by the 
time dependence of the total density matrix, of a sys- 
tem made by N noninteracting qudits each under the 
influence of a local noise when the single qudit dynam- 
ics is known. This procedure is a generalization of the 
one developed previously for a system made by only two 
parts [24J and it is valid for any initial condition of the 
total system. Moreover, it results effective in simplifying 
the work typically required by the alternative approach 
where differential equations for the density matrix ele- 
ments of the total system are solved [l3 . [25} . It is also 
an useful tool in obtaining the corresponding dynamics 
for any part of the total system described by the reduced 
density matrix, as for example in obtaining the two-body 
dynamics from a three-body system. It remains of inter- 
est how to extend this procedure to the case when direct 
interaction between the systems is present [l5| . 

Our procedure has been applied to a system of two 
identical independent qubits, each of them locally inter- 
acting with a bosonic reservoir. The general expression 
of the matrix elements of the reduced two-qubit system 



has been obtained in terms of three functions of time 
whose explicit form depends on the detailed properties 
of the reservoir but is independent on the initial state 
(Sec. mil) . These results have been used to study the 
entanglement dynamics of two identical qubits in a non- 
Markovian environment at zero temperature, taking as 
initial states the mixed EWL states defined in Eq. (fT5|) 
(see Sec. lIVp. This constitutes an extension of our previ- 
ous results obtained for initial pure Bell-like states [23 |. 
We have shown that the dynamics of entanglement, de- 
scribed by concurrence, presents revivals after periods of 
complete entanglement disappearance for some range of 
parameters, typically for high degree of entanglement of 
the initial states, as already found in the case of initial 
Bell-like states. Moreover, the phenomenon of ESD may 
occur depending on the values of parameters like purity 
or degree of entanglement of the initial state. Another 
parameter which is relevant in determining the entangle- 
ment dynamics is the degree of non-Markovianity in the 
system, that is represented by the ratio between environ- 
ment effective spectral width and the spontaneous decay 
rate of the single qubit (Sec.|V|. The revivals of entan- 
glement obtained for our two-qubit system differ from 
the ones previously obtained in the presence of interac- 
tion among qubits or because of their interaction with a 
common reservoir [isl [l6l , [TtI . The physical conditions 
examined here are similar to those typically considered in 
quantum computation and information, where qubits are 
independent and may interact with non-Markovian envi- 
ronments typical of solid state micro devices |47|]. Our 
results show that also when initial states are not pure the 
occurrence of entanglement revivals can extend the time 
when entanglement can be usefully exploited. 

We have also applied our procedure to the two-qubit 
system in a Markovian environment at finite temper- 
ature (Sec. EH, again starting from extended Werner- 
like states, generalizing thus the results present in liter- 
ature We have shown that, in the low temperature 
case (fcsT <^ hujo), ESD always occurs, at variance with 
the case of zero temperature where for some pure initial 
states it does not occur [l^ . For qubits implemented 
by quantum dots the condition may occur that, instead, 
it is applicable the high temperature limit {ksT ^ fnoo). 
Under this condition ESD always occurs and we have 
obtained a simple analytical expression of the ESD time 
in terms of the physical parameters such as purity and 
degree of entanglement of the initial state, that is a gener- 
alization of the results previously obtained for initial Bell 
states In particular we have shown that, for some 

values of the temperature, the ESD time can be shorter 
than the single qubit decoherence time and this can give 
a limitation on the time of usefulness of entanglement in 
devices made by quantum dots. 
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